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A novel interplay between d-wave superconducting order parameter symmetry and the underlying 
Cu 3d orbital based electronic structure of Copper-Oxides leads to a striking anisotropy in the 
superfluid response of these systems. In clean tetragonal materials the c-axis penetration depth 
increases as T 5 at low temperature T, in contrast to linear T behavior in the afe-plane. Disorder is a 
relevant perturbation which causes all components of the superfluid response to depend quadratically 
on temperature at low temperature. However, the cross-over temperature scale from the intrinsic 
d-wave behavior to the disorder dominated behavior for the in-plane response may be different from 
that for the out-plane response. 



PACS numbers: 74.25.Nf,74.20.Mn 

A variety of experimental data on Copper-Oxide Su- 
perconductors are compatible with a pairing state of 
d x 2_ y 2 symmetry in Cu02 planes. The existence of en- 
ergy gap nodes in this state has profound consequences 
for the low temperature electro-magnetic response. The 
in plane penetration depth increases linearly with T at 
low T in contrast to the activated behavior of conven- 
tional superconductors. The linear-T result holds in the 
regime where the response is governed by coherently 
propagating fermion quasi-particle excitations, i.e. at 
temperatures r C T « T c where Y is the superconduct- 
ing quasiparticle scattering rate. Linear-T behavior in 
the a6-plane has now been established in both chain and 
non-chain Copper-Oxides fj]-f§. Experiment 
however, indicates a much weaker T dependence of c- 
axis penetration depth. Straightforward extension of the 
d-wave model to 3D leads to line gap nodes on a warped 
cylindrical Fermi surface. This gives a linear-T c-axis 
response in apparent contradiction to experiment 0. 

In this paper we propose that weak T-dependence of 
the c-axis superfluid density p° s is a general feature of non- 
chain cuprates with crystal tetragonal symmetry, whose 
origin lies in the simple but unusual electronic structure 
of these materials. Our analysis is based on the assump- 
tion that coherent Bloch bands along c-axis are present 
at low temperatures, i.e. k z is a good quantum number. 
As discussed below, an essential feature of the electronic 
structure is that the c-axis hopping integral fj_(fc||) is a 
function of the in plane momentum fcii and is vanish- 
ingly small when fcii lies along the zone diagonals of the 
2D Brillouin zone ||||. The nodal lines of the d x i_ y i- 
wave gap A(fc||) and the zeros of t±(k\\) therefore coin- 
cide. This weakens the d-wave node contribution along 
the c-axis relative to the a6-plane and leads to a weak 
T-dependence of the c-axis superfluid response p c s . At 
sufficiently low temperatures disorder effects are always 
relevant, and tend to enhance the T dependence of p c s . 

Let us first consider the electron structure in high- 
T c compounds with tetragonal symmetry. For simplic- 
ity we consider a tetragonal monolayer system, for ex- 
ample HgBa2Cu04+5 EfJ. LDA band structure calcu- 



lations reveal that a minimal tight-binding model for 
high-T c compounds involves nearest-neighbour hopping 
between Cu 3d x 2_ y 2 and 4s and O p x and p y orbitals ||,[! . 
Thus a model Hamiltonian for HgBa2Cu04+a, expressed 
in terms of bonding, a k = 2(cosfc x p x + cosk y p y )/ojk 

where uj k = 2 J cos k x + cos k y , and non-bonding, (3k — 

2(cosk y p x — cosk x ~p y )/ujk, representation of orthogonal 
oxygen Wannier orbitals with respect to the Cu d x 2_ y 2 
orbital , is given by 
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where C k = (d k , s k , a k , flk), {Ed, E 8 ) are the energy 
of (3d, 4s) orbital with respect to the O orbital, and 
{t P d, tps, t ss ) are wavefunction overlap integrals among 
d, s, and p orbitals. The symmetry functions appear- 
ing here are fi k = 2(cosfc a: — cos ky)/ujk and z k = 
Scosk x cos k y /ujk- The Hamiltonian ([[]) differs from that 
in Ref. Q purely by retention of a strong Coulomb in- 
teraction U in the Cu d orbital. This produces a charge 
transfer insulator at half filling, a feature absent in LDA 
calculations. 

In the strong coupling limit, the four band model (|l|) 
can be simplified as a one-band model by eliminating the 
highlying Cu 4s orbitals and the highlying d — p spin 
triplets. This projection procedure can be carried out 
in two steps: firstly, by elimination of the highlying 4s 
orbital; and secondly, by solving the correlation prob- 
lem within the unit cell and treating intercell hop ping a s 
a degeneracy lifting perturbation following Refs. |12| , [13[ , 
This yields the well-known t-J model with both intra- 
and interlayer hoppings and exchanges. We assume the 
non-bonding orbitals are sufficiently lowlying that they 
are filled and inert. We notice that the dominant c-axis 
overlap is via large radius s orbitals Thus terms to 
second order in E7 1 must be retained to describe effec- 
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tive c-axis hopping via O bonding orbitals. The lowest 
energy states within the CuOi unit cell are spin 1/2 dou- 
blets | lcr) with one hole per unit cell and Zhang-Rice 
singlets |2) with two holes per unit cell. The a6-plane 
hopping matrix elements for Zhang-Rice singlets involve 
both d-a and effective a-a overlap. However along the 
c-axis, the only channel for hopping is via the bonding 
O orbitals. It can be shown that the projected c-axis 
hopping integral has the form t±(k\\) = t 
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where 



t° ± = t ss (t pd t ps M /E d E s ) 2 and M = (lo-|a+|2> is the 
overlap of the state produced by eliminating the hole 
in the O bonding orbital and the ground state doublet. 
Thus t±(k\\) contains a fen dependent factor p 2 , which 
vanishes when k x — ±k y |l4]. This is a robust fea- 
ture of realistic strong correlation models of non-chain 
tetragonal materials and directly reflects the matrix el- 
ement between Cu 4s and bonding O a orbitals. It is 
worth pointing out that afr-plane anisotropy observed in 
c-axis magneto-transport experiments may be dominated 
by precisely this effect fill . 

Now let us consider the effect of this fc|| dependent t± 
on the c-axis superfiuid response. For a clean system 
the electromagnetic response is governed by coherent 3D 
fermion quasiparticle excitations at low temperature. In 
this case the superfiuid density (or inverse square of pen- 
etration depth) is given by 



E 



k L 



dk L . 



df(X k ) _ d 2 e k e k 
d\ k dk 2 A fe 



tanh(/3A fe /2) 



(2) 



where e k is the energy dispersion of electrons, — 
yje'f, + A' 2 ,, and f(\ k ) is the Fermi function. When 
t±(k\\) — t^p 2 . it is straightforward to show that at low 
temperature the c-axis superfiuid response is 
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pt(T) ~ i iV(0)(t ( [) 2 

+ o(T 5 )) , (3) 

where iV(0) is the normal density of states on CuO 
planes, A(T) is the energy gap, A = A(0), and £1 is 
the energy cut off which is much larger than A. The 
T 5 term is from the first term in (|^): one T is from the 
contribution of the linear density of states of the d-wave 
state and the other T 4 is due to the (cosfca, — cosk y ) 4 
factor in The second term in (||) comes from 

the second term in (^), at low temperature it con- 
tributes a positive T 3 term to p c s since the d-wave gap 
A(T)/A ~ 1 - (T/A ) 3 . However, the coefficient of this 
T 3 term is proportional to (Ao/r2) 2 , which is extremely 
small. Therefore it is likely that the dominant tempera- 
ture dependence of p% should be T 5 at low temperature. 
This T 5 behavior in p c s is clearly very different from the 
— constant case where p c s varies linearly with T 
at low temperature. In Fig. EL we compare the tempera- 
ture dependence of p c s (T) / p c s (Q) in the whole temperature 



regime for the above two cases. The temperature depen- 
dence of Ps(T)/ p s (0) along a-axis is also shown in Fig. 
(|l|) for comparison. 




FIG. 1. Normalized superfiuid den- 

sity p s (T)/ p a (0) vs T/T c . The energy dispersion of electrons 
e k = — 2t(cos k x + cos k y ) — At' cos k x cos k y — 2t±(k\\ ) cos k z 
is used in the calculation. (a) is for the model with 
= i±(cosfca; — cosfc H ) 2 /4 along c-axis. (b) and (c) are 
for the model with t±(k\\) = ij_ along c- and a-axis, respec- 
tively. t = l,t' = -0.25, t° ± = 0.1, T c = 0.14, and the filling 
factor is 0.4. (The results are not sensitive to the change of 
parameters) . 

There are several effects which may lead to a finite hop- 
ping integral along the c-axis in the vicinity of the gap 
nodes in clean systems and thus a stronger T dependence 
of p c s . Many non-chain materials have some degree of or- 
thorhombicity, i.e. t x ^ t y . Clearly this will affect the 
location of the gap nodes and the zeros of t± (k\\)- Follow- 
ing the previous treatment for the model ([!]), it can be 
shown that the projected c-axis hopping integral is still 
a function of fen, but tj_(k\\) now has approximately the 

form tj^cosfc^ — -r- cos k y ) 2 /A. Thus the zeros of t±(k\\) 
are no longer on k x = ±k y planes. The gap nodes may 
not lie on k x — ±k y planes either, but there is no symme- 
try reason requiring the gap nodes to coincide to the zeros 
of tj_(ki\). Thus t± is now finite at the gap node. Direct 
interplane 0-0 overlap can also lead to a finite dispersion 
around the gap nodes and eliminate the zeros in t±(k\\) 
altogether. In either case, we may model the c-axis elec- 
tronic structure by = t\n\ + t n ° de with t n ° de the 
value of t± at the gap nodes. Both effects are expected 
to be weak, but they generate a small linear term in p c s , 
with a slope proportional to (t™ de ) 2 , which dominates 
the c-axis superfiuid response when T « t™ de . 

Sufficiently close to a nodal line the c-axis hopping rate 
always becomes small in comparison with the impurity 
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scattering rate, and thus disorder is always a relevant 
perturbation. Disorder introduces a finite quasiparticle 
lifetime and consequently a finite density of states at the 
Fermi level. At low temperature the correction to the su- 
pcrfluid density caused by finite quasiparticle scattering 
rates is approximately given by 
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It can be shown that both the ab plane and the c-axis su- 
perfluid responses behave as T 2 when T is much smaller 
than the quasiparticle scattering rate T jL7|. However, 
the coefficients of the T 2 terms are very different for 
the ab plane and the c-axis responses. On the ab plane, 
pf (T)/pf (0) ~ -(A /r)(T/A ) 2 ; while along the c- 

axis, p c s {T)/p c s {0) (r/A )(T/A ) 2 . Thus the T 2 

behavior of the c-axis response is by a factor (T/Aq) 2 
weaker than that on the ab plane (generally r < Ao) 
due to the p\ factor in t_i_(fe|i). Notice that disorder 
has opposite effects in ab and c-axis; it weakens the In- 
dependence of the in-plane response and strengthens the 
T-dependence of the out of plane response. 

Another important aspect of disorder is that it may 
disrupt the symmetry of bonding orbitals about the Cop- 
per site. The selection rule preventing hopping for mo- 
mentum along k x = ±k y is then relaxed. For example, 
consider the effect of interlayer defects on O site. This 
leads to a random component of t± due to random vari- 
ation in the matrix element Mq for a Zhang-Rice sin- 
glet discussed earlier. Near the gap nodes this fluctuat- 
ing component in tj_ dominates c-axis hopping, and thus 
electrons in the vicinity of nodes are well described by an 
impurity assisted hopping model which has been intro- 
duced in the literature on phenomenological grounds . 
Impurity assisted hopping gives a new conduction chan- 
nel and has a direct contribution to the c-axis superfluid 
density p c s<imp : 
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Thus the superfluid response induced by the impurity as- 
sisted hopping p c s imp depends on the impurity average of 
the assisted hopping matrix element (St±(q)St±(—q)) imp . 
If the impurity scattering is isotropic in space, i.e. 
(6t±(q)5t±(—q))i mp is independent on q, then the as- 
sisted hopping contribution to p c s in a d-wave state is zero, 
resulted from the vanishing average of the d-wave gap on 
Fermi surface. If, however, the scattering is anisotropic, 
for example if (St±(q)St±(—q))i mp has the Lorentzian 



strong forward scattering limit (Sk 
as T 2 at low temperature: 
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where Tj_ = 2ttvqN(0). For a finite Sk, p c s imp behaves 
similarly to the Sk — > case, but the overall amplitude 
of p c s ^ mp decreases with increasing Sk. 

The above discussion indicates that both the lifetime 
effect and the impurity assistant hopping contribute a 
T 2 term to p c s at low temperature. The T 5 behavior 
of p c s from the coherent tunneling of electrons is there- 
fore observable only when T c ^> T ^> max(r,T*), where 
T* ~ ip^Ao/^) 2 ] 1 / 3 !^ is a characteristic temperature 
above which the contribution from the coherent hopping 
of electrons to p c s is larger than that from the impurity 
assisted hopping. For a highly anisotropic and clean high- 
T c superconductor, generally T± <C f j_ and t° ± is smaller 
than or of the same order as Ao. However, because of 
the 1/3 exponent in T*, T* is a rather large temperature 
scale compared with even T. Thus assisted hopping may 
affect the c-axis superfluid response in a larger tempera- 
ture region than the scattering lifetime effect. 

Summarising the discussion above, we may draw a 
qualitative phase diagram for the temperature depen- 
dence of p% and pf at low temperature: When T <C T, 
the lifetime effect is important, both p ab and p c s vary 
quadratically with T. When T* > T > T, the im- 
purity assisted hopping becomes more important to the 
c-axis superfluid response. In this case, p c s varies still 
quadratically with T, but the coefficient of the T 2 term 
in p c s (T) I Ps(0) is proportional to r_i_/(i^) 2 A, which is 
generally smaller than that in the T T case where the 
coefficient of T 2 term in p c s (T) / p c s (Q) is proportional to 
T/A 3 . When T c > T > T for pf or T c > T > T* 
for p c s , the disorder effect is not so important and the 
temperature dependences of p c s and pf in a clean sys- 
tem should be recovered. Thus there is a crossover from 
T 2 to T when T - T for pf and from T 2 to T 5 when 
T ~ T* for p% as T increases at low temperature @. 
A simple formula which gives a qualitative measure for 
the crossover of p c s from the intrinsic T 5 behavior to the 
disorder T 2 behavior is p c s (T) m p c s (0) + aT 2 [l + {T/T*) 3 } 
with a a scattering potential dependent constant. The 
above picture may be also valid for the case of substitu- 
tional impurities such as Zn, but we cannot rule out a 
possible formation of a "impurity band" which may have 
a stronger effect on the low temperature behavior of p c g 
than the disorder effects previously discussed. 

In comparing the predictions of this model to experi- 
ment for a given system, it is essential to assess whether 
coherent Bloch band are present at all. The answer to 
this question in the normal state of highly anisotropic 
copper-oxides, such as BSSCO, appears to be no. The 
c-axis mean free path is generally very short compared 
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with the c-axis lattice constant. Microwave |19[| and ther- 
mal Hall measurements J2(]], however, reveal that the 
mean free path grows very rapidly in the superconduct- 
ing state. The superconducting quasiparticle scattering 
rate is therefore much lower than the extrapolated nor- 
mal state scattering rate. Thus it is our belief that the 
available experimental data favor coherent Bloch band 
formation along c-axis at low temperature below T c . For 
BSSCO, tj_ is small and p c s might be dominated by the 
impurity assisted hopping [Q. However, for other com- 
pounds with lattice tetragonal symmetry and relatively 
lower anisotropy, such as LaSrCuO, TIBaCaCuO and Hg- 
BaCuO, we believe that the coherent band should have 
substantial contribution to p c s . Experimentally, a cru- 
cial test for this can be done by carefully measuring and 
analysing p c s data in these compounds to see if the pre- 
dicted T 5 behavior exists. (The T 5 behavior in p c a is not 
applicable to YBCO as it has CuO chain layers and non- 
fiat CuO planes fl|.) 

In conclusion, we have demonstrated that there ex- 
ists a remarkable interplay between underlying electronic 
structure and orbital parameter symmetry in Copper- 
Oxides. The c-axis electromagnetic response reflects not 
only the symmetry of the order parameter but also the 
relative symmetry of the Copper orbitals involved in c- 
axis and afe-plane conduction. At low temperature the 
disorder effect dominates and both p c s and p^ b have T 2 
dependence. Above some characteristic temperatures, 
which may be different for p c s and p° b , the temperature 
dependences in a clean system, linear T for p^ b and T 5 
for p c 8 , are expected. 

The authors are grateful for valuable discussions with 
D. C. Morgan, N. E. Hussey, and J. R. Cooper of the 
IRC in Superconductivity at Cambridge. 
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